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Late growth stages and post-growth diffusion in organic epitaxy: PTCDA on Ag(lll) 

B. KrauseQ A. C. Diirr, F. SchreiberQ and H. Dosch 
Max-Planck-Institut fur Metallforschung, Heisenbergstr.3, 
70569 Stuttgart and Institut fur Theoretische und Angewandte Physik, 
Universitat Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany 



o 
o 

(N 

Oh' 



o 



c 

o 
o 



> 

(N 
O 

o> 
o 
^- 
o 

c3 



c 

o 
o 



13 



0. H. Seeck 

IFF, Forschungszentrum Jiilich GmbH, D-52425 Julich, Germany 
(Dated: February 2, 2008) 

The late growth stages and the post-growth diffusion of crystalline organic thin films have been 
investigated for 3,4,9, 10-perylenetetracarboxylic dianhydride (PTCDA) on Ag(lll), a model system 
in organic epitaxy. In situ x-ray measurements at the anti-Bragg point during the growth show 
intensity oscillations followed by a time-independent intensity which is independent of the growth 
temperature. At T > 350 K, the intensity increases after growth up to a temperature-dependent 
saturation value due to a post-growth diffusion process. The time-independent intensity and the 
subsequent intensity recovery have been reproduced by models based on the morphology change 
as a function of the growth temperature. The morphology found after the post-growth diffusion 
processes has been studied by specular rod measurements. 

PACS numbers: 68.55.Ac, 81.10.Aj, 61.66.Hq 



I. INTRODUCTION 

Crystalline small-molecule organic semiconductors 
such as perylene and phthalocyanine derivatives show in- 
teresting electrical and optical properties, which can be 
exploited in electronic and optoelectronic devices based 
on thin films deposited by organic molecular beam epi- 
taxy (OMBE). However, the intrinsic properties of the 
organic crystals are closely related to the molecular and 
crystalline anisotropy, and the crystalline quality of the 
organic thin film has a strong impact on the resulting 
device performance. In order to control the properties of 
organic thin films, the understanding of the basic growth 
processes is necessary. 

While the structure and epitaxy has been thoroughly 
studied for several model systems, such as PTCDA 
QESSIiSHIi, on Ag(lll) the growth dynamics 
and evolution of the resulting morphology have been less 
systematically investigated. Importantly, even for a sys- 
tem with well-defined monolayer epitaxy, the later stages 
of the growth (after a certain threshold thickness) can, 
of course, exhibit islanding (Stranski-Krastanov (SK) 
growth) and a very rough resulting morphology. It was 
recently found that PTCDA on Ag(lll), one of the struc- 
turally best characterized OMBE systems, exhibits a rel- 
atively smooth morphology (kinetically limited) at low 
temperatures (T < 350 K), but that the growth behav- 
ior undergoes a transition to a very clear SK scenario at 
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high temperatures (T > 350 K). In this paper, we present 
a real-time study of the later stages of the growth (is- 
landing process) and the post-growth diffusion processes, 
which have a profound impact on the resulting morphol- 
ogy. We attempt a comprehensive discussion of the ef- 
fects related to the temperature dependent transition of 
the morphology. 

PTCDA (C 24 6 H 8 , see Fig. □ (a)) is an archetypal 
system for OMBE Q. It is a flat molecule which crystal- 
lizes in two monoclinic bulk polymorphs, the a and the 
j3 phase [9|, |lfj, 111], [l2| • In the following, we will always 
refer to the a bulk notation. For both polymorphs, the 2 
molecules of the unit cell order in a herringbone structure 
with the molecular plane parallel to the PTCDA(102) 
plane. At room temperature, the stacking distance be- 
tween subsequent molecular sheets is d F — 3.22 A and 
d F = 3.25 A in PTCDA(102) direction. 

On Agflll), PTCDA forms well-ordered epitaxial 
films 0, Ej ■ It has been found that the binding between 
the commensurate first monolayer and the substrate is 
much stronger than the binding between the subsequent 
layers 01 • The structure and the morphology of PTCDA 
films vary significantly with the growth temperature, T 
0,113. For T < 350 K, the films are strained and have 
a relatively homogenous thickness. For T > 350 K, the 
films consist of relaxed islands on top of two wetting lay- 
ers. For a variety of growth conditions, a coexistence of 
both bulk polymorphs has been observed. 

A detailed study of the epitaxy and the post-growth 
morphology of the PTCDA films studied here has been 
published elsewhere 0, fl5l 1 1(1 Il7| . In this paper, we 
focus on the results of specular intensity measurements 
(with the momentum transfer q = q z e z perpendicular to 
the sample surface) during and directly after the growth 
of thin PTCDA films. This type of x-ray measurement 
is very suitable for studying the evolution of the out-of- 
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plane structure and morphology of organic thin films in 
real time. It contains the laterally averaged information 
about the thin film structure in z-direction. 

Two types of experiments arc reported in the following: 
measurements of a section of the specular rod, i.e. inten- 
sity measurements as a function of q z , after the growth, 
and time-dependent intensity measurements at one spe- 
cific point of the specular rod during and after the de- 
position of the PTCDA films. For this, the anti-Bragg 
point of the film with q* = tt/cIf has been chosen (see 
below) . 

The first type of experiment gives insight into the crys- 
talline structure of the films. It will be shown how the 
growth temperature dependence of the morphology influ- 
ences the specular intensity. The second type of exper- 
iment gives time-dependent information about the ma- 
terial distribution between the different layers. It con- 
tains information about the growth process and the post- 
growth diffusion. The initial stages of the growth are de- 
scribed in detail elsewhere 0|. Here, the later stages of 
the growth (after the growth of the wetting layer) and 
the post-growth diffusion are discussed. 

The paper is organized as follows. In Sec. |H] the de- 
tails about the sample preparation are given. The exper- 
iments are described in Sec. IIIII In Sec. IIVI the results 
are discussed, and theoretical models for the observations 
are developed. The main results and conclusions of the 
paper are summarized in Sec. E] 
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FIG. 1: (a) 3,4,9, 10-perylenetetracarboxylic dianhydride 
(PTCDA) . (b) Measured (circles) and calculated (lines) spec- 
ular rod intensity as a function of the growth temperature T 
for samples with F — 0.8 — 2 A/min and d = 50 A. 



II. EXPERIMENTAL 

The organic thin films have been deposited by molecu- 
lar beam epitaxy in a portable UHV chamber (base pres- 
sure 10 -10 mbar) equipped for x-ray experiments |l8| . A 
Ag(lll) single crystal with an out-of-plane mosaicity of 
0.13 ±0.02° has been used as a substrate. Prior to depo- 
sition, the substrate has been cleaned by repeated cycles 
of sputtering (30 min at room temperature, Ar + energy 
500-600 eV, Ar pressure k 6x 10~ 5 mbar) and annealing 
at « 700 — 900 K for at least 15 min. The purified and 
thoroughly outgassed PTCDA has been evaporated from 
a Knudsen cell kept at approximately 575 K. 

Two sets of samples with the average film thickness 
d = 50 A and growth temperatures T in the range of 
197 - 473 K have been studied, one with a flux (growth 
rate) F = 0.8 - 2.0 A/min (low F), the other with F = 
12- 14 A/min (high F). 

The experiments have been performed at the Ham- 
burger Synchrotronstrahlungslabor (HASYLAB) at the 
beamline Wl with the wavelength A = 1.237 A. 



III. RESULTS 
A. Specular rod measurements 

Fig. H shows examples of the specular rod measure- 
ments including the PTCDA(102) reflection for samples 
grown at low F. The measurements have been performed 
at the respective growth temperature of the samples. In- 
dependent of T, the PTCDA(102) plane is parallel to 
Ag(lll). For T < 358 K, 1-3 Laue oscillations around 
PTCDA(102) have been observed. In this temperature 
range, the width of the PTCDA(102) peak corresponds 
to the average film thickness and is independent of T. 
For T > 358 K, the peak width of the Bragg peak de- 
creases with increasing T, and the Laue oscillations re- 
duce to two shoulders on the Bragg peak or disappear 
completely. With increasing T, the rod develops a pro- 
nounced minimum at g, » 1.4 A -1 . While the rod mea- 
surements show a strong temperature-dependence, the 
rocking scans through PTCDA(102) are independent of 
T and show the same peak width of 0.13 ± 0.02° corre- 
sponding to the mosaicity of the substrate. A detailed 
discussion of the modelling of the data is presented in 
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FIG. 2: Simulation of the specular rod of a thin PTCDA 
film on Ag(lll). At the anti-Bragg point of the PTCDA film 
(q* — 77 /(If), the scattering of subsequent layers interferes 
destructively (see inset). For layer-by-layer growth, the in- 
tensity at the anti-Bragg point oscillates as a function of the 
deposition time. 



Sec. EH 

The measurements for the high F samples (not shown 
here) show a similar behavior, but the changes in the rod 
measurements and the peak shape are shifted to higher 
temperatures by about 15 K. The Laue oscillations and 
the peak width expected for the average film thickness 
are still observed for T < 358 K, where the oscillations 
are already damped out for the low F samples. 



B. In situ experiments 

In order to follow the growth of thin PTCDA films, the 
x-ray diffraction intensity from the molecular layers has 
been measured in real-time. In kinematic approximation 
the specular x-ray scattering is the sum of the scattering 
contributions from the film and the substrate, 
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Here, Aq z is the difference between q z and the nearest 
specular Bragg peak of the substrate with q z = I2n/ds- 
dp and ds are the lattice spacings of the film and the sub- 
strate, fp and fg are the corresponding form factors, d s i 
is the lattice spacing of the topmost layer of the substrate, 
do is the distance between the substrate and the first layer 



of the film, 9 n is the fractional coverage of the fcth layer 
within the film, and erg is the Gaussian roughness of the 
substrate. At the anti-Bragg point of the PTCDA film 
(q* z = tt /dp) the contribution of the substrate (the sec- 
ond term in Eq.^l equals fp l)( fe-1 )0fc(t), which can 
be written as JfAO, where A6 is defined as 



ven (^) • 

.(2) 

This is the difference between the material (in fractions 
of monolayers) deposited in odd layers (i.e. the coverage 
of layer 1 + the coverage of layer 3 + ...) and the material 
deposited in even layers. This simplification is possible 
since at the anti-Bragg point the scattering amplitude of 
all even and all odd layers has the same phase. It has 
been assumed here that the lattice parameter of the film 
is constant across the entire film thickness. This assump- 
tion is supported by the specular rod measurements. 

The scattered intensity at the anti-Bragg point, nor- 
malized to the scattering of the substrate, is 



I(q* z ,t) = l + AAe + B{ABf 



(3) 



with 
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(4) 



The parameters A and B depend on the materials in- 
volved. For PTCDA on Ag(lll), inserting the atomic 
form factors and the lattice parameters determined from 
the rod measurements (see Sec. lIII A"fr leads to the values 
A = -0.73 and B = 0.06. 

In the case of layer-by layer growth, characteris- 
tic intensity oscillations with the minimum intensity 
Imin — 1 + A + B and the maximum intensity I max — 1 
are observed as a function of the deposition time 
EHIHE! (see Fig. 

Fig. |3| shows a typical intensity measurement during 
the deposition of a PTCDA film, normalized to the in- 
tensity scattered by the substrate measured before de- 
position. The time is given in units of the monolayer 
deposition time, tml, of the first few layers as derived 
by the observed intensity oscillations. The temperature- 
dependence of the intensity oscillations during the de- 
position of the initial monolayers has been described in 
detail elsewhere ^5|. In the following we briefly sum- 
marize the most important observations, and then focus 
on the analysis of the later stages of the growth and the 
post-growth changes. 

Before opening the shutter of the evaporation cell, the 
constant intensity, Is, scattered by the substrate is mea- 
sured. During deposition, the intensity oscillates between 
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FIG. 3: Typical intensity measurement at the anti-Bragg 
point during the deposition of a film with T = 393 K and 
F Kil A/min. 



the maximum value I(q%, 2tml) % 1 ±0.02 and the mini- 
mum value I(q*,TML) > 0.33±0.02 for t < 3tml- tml is 
the deposition time for 1 ML. For T < 358 K the intensity 
oscillations for t < 3tml are damped, for T > 358 K they 
are undamped and show no temperature dependence. In- 
terestingly, during further deposition, the intensity ap- 
proximates the asymptotic value I(q£, oo) = 0.65 ± 0.02, 
independent of T. 

After the shutter has been closed, the intensity in- 
creases with t up to a saturation value, Is a t > HqIi °°); 
indicative of post-growth structural changes. Fig. 0] (a) 
shows this intensity recovery as a function of the growth 
temperature. For T < 358 K, I(q*,t) = I(q*,oo) = 
0.65 ± 0.02 remains constant after deposition, while for 
T > 358 K, the intensity increases as a function of the 
time after interrupting the deposition. The saturation 
intensity shown in Fig. 0] (b) increases with T up to 
Isat = 1 ± 0.03 , while the characteristic time, Ts a t, for 
the saturation process decreases with T (see Sec. II VI for 
a quantitative analysis). 



IV. DISCUSSION 



For the PTCDA films discussed here, a morphology 
change from a smooth or facetted film to islands on a 
wetting layer has been observed as a function of T |2Hl5|. 
The transition to the island morphology has been found 
around T « 350 K which is comparable to the tempera- 
ture where the changes in the rod measurements and in 
the in situ intensity measurements have been observed. 
This suggests that the observations reported here are re- 
lated to the morphology change. 
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FIG. 4: (a) Intensity at the anti-Bragg point after the deposi- 
tion of a film as a function of T, measured for several samples 
with d at 50 A and F = 0.8 - 2 A/min. The shutter of the 
evaporation cell has been closed at t = 0. (b) Saturation 
intensity as a function of T. 



A. Rod measurements 

Qualitatively, the observed features of the rod mea- 
surements correspond to the following scenario. The 
PTCDA films grown at low temperatures (for the low 
F samples, this corresponds to T < 358 K) are rela- 
tively smooth as indicated by the Laue oscillations. The 
films cover essentially the entire substrate surface since 
the width of the Bragg peak corresponds to the aver- 
age film thickness. For higher temperatures, the Laue 
oscillations vanish indicating an increasing roughness of 
the film. The decreasing width of the Bragg peak cor- 
responds to high islands, and the intensity minimum at 
q z w 1.4 A -1 indicates a 2 ML thick wetting layer be- 
tween the islands. The observations agree with the mor- 
phology observed by AFM 0, implying that the entire 
PTCDA film is ordered. 

For all films, only a single PTCDA out-of-plane Bragg 
peak around q z ~ 1.95 A -1 has been observed, even if the 
in-plane measurements have confirmed the coexistence of 
the a and the (3 phase [2|. The expected splitting be- 
tween the (102) Bragg reflections of both polymophs is 
AQ = 2ir/df - 2tt/<4 « 0.02 A" 1 . For T < 358 K, 
AQ is much smaller than the observed width of the 
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FIG. 5: Schematic of the morphology assumed for the simu- 
lations, (a) Low temperature model for T < 358 K, (b) high 
temperature model for T > 358 K. 



PTCDA(102) reflection and cannot be resolved. For the 
narrow peaks observed at higher T, the peak splitting 
is not observed due to the dominance of the a phase at 
these temperatures. 

Interestingly, the narrowing of the PTCDA(102) peak 
as a function of the growth temperature reported here has 
also been observed for PTCDA grown on H-passivated 
Si(100) H2, for PTCDA on Au(lll) @, and for annealed 
films of PTCDA 0. 



B. Modelling of the rod measurements 

Based on the morphological features observed with the 
AFM [3 , the following model for the X-ray data has been 
developed. 

In kinematic approximation the specular x-ray scatter- 
ing intensity of each mosaic block of the substrate cov- 
ered by a thin film is I(q z ) = \(S(q z )}\ 2 with S(q z ) as 
described in Eq. ^ Since the morphology discussed in 
this section is observed long after the growth process and 
therefore time-independent, the parameter t of Eq. ^has 
been omitted. 

Fig. [5] (a) shows a schematic of the model applied to 
the rod measurements with a temperature-independent 
peak width (comparatively smooth film). Two types of 
mosaic blocks have been introduced to describe this low- 
T morphology: (1) mosaic blocks with np film layers and 
a Gaussian roughness <7p, described by 



= 0.5 



1 - erf 



np 



OF 



(5) 



and (2) mosaic blocks covered by a film of only a few 
monolayers and 



> for k < 3 
for k > 3 



(G) 



where k is the layer index. The total intensity scattered 
from the illuminated area is 



I(q z ) = (l-A 2 )\(Si(q z ))\ 2 +A 2 \(S 2 (q z ))\< 



(7) 



where Si(q z ) and S 2 (q z ) are the scattering amplitudes of 
the different types of mosaic blocks, and A 2 is the fraction 
of the substrate surface covered by the type- (2) mosaic 
blocks. 

For the rod measurements with a decreasing peak 
width (islands on a wetting layer, high-T films), the 
model is slightly modified (see Fig. [SI (b)). The islands 
(i.e. the PTCDA film of the type (1) mosaic blocks) are 
assumed to be molecularly flat with 



for k < n 
for k > n. 



(8) 



They have a Gaussian height distribution with the av- 
erage height np and the standard deviation a p. The 
scattered intensity of the film is 
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(9) 



where A 2 is the area covered by the type- (2) mosaic 
blocks corresponding to the wetting layer, and S\(q z ,n) 
and S 2 {q z ) are the scattering amplitudes of the different 
types of mosaic blocks. 

These two models describe all specular X-ray mea- 
surements. The results of the calculations are shown as 
straight lines in Fig. ^ Fig. [S] summarizes the parame- 
ters used for the simulations. For T > 358 K, the lattice 
parameter varies as expected from the thermal expan- 
sion of PTCDA in (102) direction PJ]. For T < 358 K, 
the effective lattice parameter seems to be approximately 
temperature independent instead of decreasing further. 
A possible explanation for this is that the ratio of the 
PTCDA (3 phase with a slightly larger lattice parameter 
than the a phase increases for T < 358 K as already in- 
dicated by in-plane measurements of thin PTCDA films 

a- 

The area covered by the wetting layer and the island 
height increase with T for T > 358 K, as observed in 
a similar AFM study 0- Since the island density in- 
creases with increasing deposition rate and decreasing 
growth temperature, the islands grown at high F are al- 
ways smaller than the islands grown at the same T but 
lower F. 
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FIG. 6: Parameters used for the simulations shown in Fig.Q 
Open symbols: low F data, filled symbols high F data, (a) 
Lattice parameters of the PTCDA film as a function of T, (b) 
area covered by the type-(2) mosaic blocks corresponding to 
the wetting layer, and (c) height uf of the type-(l) mosaic 
blocks corresponding to areas covered by a thicker PTCDA 
film. The error bars of uf indicate of, i.e. the roughness in 
the low temperature model and the island height distribution 
in the high temperature model, respectively. 



The thickness of the wetting layer is approximately 
2 ML and the lattice spacing between the topmost layer 
of the substrate and the first layer of the PTCDA film is 
d a = 2.8 ± 0.1 A, independent of T. 

For PTCDA on Au(lll), similar measurements of the 
specular rod have been performed [8j. These measure- 
ments have been modeled as the coherent scattering of 
a film with a laterally averaged density. In contrast to 
this, here the incoherent addition of the intensity has 
been chosen since the AFM images of PTCDA/Ag(lll) 
show that the morphology varies on a length scale much 
larger than the average mosaic block size of 500 A, and 
the scattering of different mosaic blocks adds up incoher- 
ently. 



C. In situ growth experiments 

In the following, two aspects of the real-time growth 
observations will be discussed in detail, the asymptotic 
intensity, I(q*, oo), and the intensity increase observed 
immediately after the deposition. Both features are re- 
lated to the morphology of the films. 
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FIG. 7: Schematic explaining the asymptotic intensity. The 
roughness of the surface is assumed to be Gaussian. The 
average sample surface corresponds to the layer mo . Partially 
filled layers k in the same gray scale add up to a competely 
filled monolayer. 



1. Asymptotic intensity 

For all samples, a constant asymptotic intensity 
J(g*,oo) has been observed after the initial growth 
stages. With Eq. the experimentally observed value 
I(q*, oo) = 0.65 ± 0.02 corresponds to A8 = 0.5 ± 0.03. 

In the following a model for this asymptotic intensity 
is developed. For simplicity only the coherent scattering 
from one mosaic block is discussed. Strictly, this model 
only applies to the case where the roughness of the film is 
much smaller than the length scale of the mosaic block, 
as observed for the later growth stages of the relatively 
homogenous films grown at T < 350 K. At higher growth 
temperatures, the scattering from the PTCDA films after 
the initial growth stages is a superposition of the inco- 
herent scattering of different types of mosaic blocks, as 
already discussed in Sec lIV Bl However, the incoherent 
addition does not change significantly the value for the 
asymptotic intensity, and the calculations can be per- 
formed in analogy to the coherent case. 

For a Gaussian height distribution (which has been 
found on a length scale of the x-ray coherence length, 
approximately 1 ^m, for films with T < 350 K after 
growth Q), the coverages Ok of a film with an average 
layer number mo are determined by (see Eq. 



h= 2{ l - eTf 



(10) 
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(Fig. UJ. For sufficiently large of, Eq. ^]can be charac- 
terized by 



'rao+A 



and 



= l-6» 
= 0.5. 



mo — A 



(ii) 



A function with these features leads to AO rj 0.5 and 
thus I(q%, oo ) = 0.65, which is shown in the Appendix. 

For the system PTCDA/Ag(lll), two cases have to be 
distinguished: The growth at T < 350 K, and the growth 
at T > 350 K. 

For T < 350 K, the layer- by-layer growth of the ini- 
tial layers (indicated by the oscillations) changes with 
increasing deposition time to a Gaussian height profile 
explaining the asymptotic intensity during and the spec- 
ular rod measurements after deposition (see Sec. IIV Bl . 
This case is consistent with the above model. 

At T > 350 K, separate islands on top of a wetting 
layer are observed after growth. This morphology de- 
viates clearly from a symmetric roughness profile. The 
dominant part of the sample is covered by two monolayers 
(AO = 0), corresponding to I(q* z ,t) = 1. The islands give 
only a negligible contribution to the measured intensity 
since they cover only a small part of the area. In addition 
to this, only the topmost layer of the island contributes 
to the scattering since the scattering of all lower layers 
interferes destructively at the anti-Bragg point. How- 
ever, one possible explanation for the observed intensity 
corresponding to AO — 0.5 is that during growth the 
third layer is partially filled with a coverage of approx- 
imately 0.5 ML. This coverage is not stable, and after 
growth, the molecules diffuse to the islands and only the 
2 wetting monolayers observed in the rod measurements 
remain. In the next section, we will see that this model is 
consistent with the observation of the intensiy recovery 
after growth. 



2. Intensity recovery after deposition 

For samples deposited at growth conditions corre- 
sponding to islands on a wetting layer, a post-growth 
intensity recovery has been observed, while for sam- 
ples grown at lower temperatures the intensity does not 
change after growth (on accessible time scales of several 
hours). During growth, the coverage changes of the dif- 
ferent layers are due to deposition and diffusion. The 
post growth intensity changes are a sign for interlayer 
diffusion after growth (desorption is a negligible effect in 
the temperature range studied here). 

The following model is consistent with the temperature 
dependence of the intensity recovery. For T < 350 K, 
after growth the samples have a Gaussian roughness on 
a scale much smaller than the mosaic block size. This 
roughness is not significantly changed by surface diffusion 
because T is too low. Therefore, the scattered intensity 
does not change after growth. For T > 350 K, post- 
growth diffusion of molecules from the top of the wetting 
layer to the islands takes place. 
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FIG. 8: (a) Area covered by the wetting layer as a function 
of T. Open symbols: fit result of the post growth intensity 
changes; filled symbols: result of the fit of the specular rod 
measurements, (b) Characteristic time for the post growth 
intensity changes as a function of T. 



To model the x-ray scattering of the inhomogeous mor- 
phology at T > 350 K, the scattered intensity of 3 
types of mosaic blocks is added incoherently: (1) Mo- 
saic blocks belonging to islands with even layer num- 
ber, the coverage difference AO = 0, and the normal- 
ized scattering intensity I even = 1 according to Eq. [3J 
(2) mosaic blocks belonging to islands with odd layer 
number, the coverage difference AO — 1, and the nor- 
malized scattering intensity I dd = 1 + A + B; (3) mosaic 
blocks covered by the wetting layer and some additional 
molecules, with the time-dependent scattering intensity 
I 3 (t) = 1 + AAO + B(A0) 2 and the coverage difference 
A0ML(t) = 0.5e _t / TSat covering the area A3. Ts a t is an 
effective parameter comprising the surface diffusion bar- 
rier, the size of the critical nucleus, and the temperature 
dependence of the surface diffusion. We assume that the 
islands with even and odd layers each cover half of the 
remaining area, which is 0.5(1 — A3). The total scattered 
intensity is 



I(t) = A 3 I 3 (t) + 0.5(1 - A 3 )(I even + I odd ). 



(12) 



Fig.EJshows the fit results of the intensity measurements 
after deposition. The area covered (only) by the wetting 
layer (Fig. El (a)) increases with the deposition tempera- 
ture up to nearly the entire sample surface. This is con- 
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sistent with the result from specular rod measurements, 
also shown in Fig. 0(a), and the AFM observations 0- 
The characteristic time, Ts a t, for the recovery process is 
shown in Fig. 0(b). It decreases with increasing T, as ex- 
pected for a diffusion-related process, from Ts a t ~ 200 s 
at T = 373 K to T Sat « 50 s at T = 473 K. 

Similar to the x-ray intensity recovery after growth re- 
ported here, an intensity recovery after the growth of 
GaAs has been observed by Reflection High Energy Elec- 
tron Scattering ^ n ^is case, the recovery cor- 
responds to a smoothening of the surface after the in- 
terruption of the Ga flux due to the post-growth surface 
diffusion. In contrast to this, the post-growth intensity 
recovery of PTCDA is explained as the result of the in- 
terplay between post-growth surface diffusion processes 
and the morphology of the PTCDA films. 



SUMMARY AND CONCLUSIONS 



The later growth stages, the post growth diffusion, and 
the structure of thin films of PTCDA/Ag(lll) have been 
studied. The specular rod data are strongly influenced 
by the temperature-dependent morphology of the films. 

The thickness of the wetting layer (2 ML), the area 
covered by the islands, and the binding distance between 
the first monolayer and the substrate have been deter- 
mined. In-situ measurements at the anti-Bragg point 
have been used to derive information about the growth 
of PTCDA. After the intensity oscillations during the 
deposition of the wetting layer (indicating layer-by-layer 
growth), a constant intensity has been observed, followed 
by an increase of intensity after growth for T > 350 K. 
The measurements have been analyzed quantitatively us- 
ing the following growth model: At T < 350 K, after the 
deposition of the wetting layer the growth mode changes 
to 3D growth with a Gaussian roughness which is stable 
after growth. At T > 350 K, the growth of the wetting 
layer is followed by island growth. During growth, the 
third layer is partially filled. After growth, due to the 
strong surface diffusion this material is also incorporated 
into the islands, and the temperature dependence of the 
time scales of the post-growth diffusion has been deter- 
mined. 

Morphology changes as a function of the growth tem- 
perature have a strong impact on the performance of op- 
toelectronic devices. The observed post-growth diffusion 
may be of special interest for multilayer structures since 
the interface quality can be improved by growth inter- 
ruptions allowing for the incorporation of residual surface 
molecules. Also, the resulting island morphology might 
be exploited similar as in the Stranski-Krastanov growth 
of inorganic semiconductors. 
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APPENDIX A: ASYMPTOTIC INTENSITY 

It is shown that a function 0k with the features de- 
scribed in Eq.^] leads to AO as 0.5 and thus I(q'*, oo) — 
0.65. AO can be written as 



AO = j^kt-l)*- 1 



k=l 

mo — 1 oo 

= o k (-i) k -'+ £ ^(-l)*- 1 

k—1 fc-mo+1 

+e mo {-i) m °-\ (Ai) 

Here, the sum has been separated into the contribution 
of the layers with k < too, the layers with k > too, and 
the contribution of the layer mo. With k = too — A, we 
get 



Atf = ]r mo _ A (-ir°- 

A=l 



A-l 



+ J2 ^„o+a(-1)" 1o+a - 1 + 9 mo (-l) m °- 1 (A2) 

A=l 

With Eq.Eland (-l)™o-A = (•_ 1 )m +A it f u ows 

mo— 1 

ao = J2 (i-^ 0+ A)(-ir o+A - 1 

A=l 
oo 

+ £ e mo+A (-i) m »+ A - 1 + mo (-ir°- x 

A=l 

mo— 1 

_ (_l)mo+A-l 

'"'"^ A " ' ^6> mo (-l) m °-(A3) 



A=l 



+ £ mo +A(-i) r 

A=m + 1 

For too ^> a pre da, and with mo = 0.5 this reduces to 

C mp— l \ 
£(-l) A + 0.5 1=0.5 (A4) 

for both mo even and to odd. 

The physical reason for this is the following: If enough 
layers contribute to the Gaussian surface roughness, a 
layer Too can always be chosen with mo ks 0.5, as indi- 
cated in Fig. [7| The X-rays scattered by the layer too + 1 
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interfere constructively with the X-rays scattered by the 
layer to — 1. 

If the coverage of the layer too + 1 is shifted to the 
layer mo — 1, the specularly scattered intensity at the 
anti-Bragg point does not change. Because of the sym- 
metry of the Gaussian lineshape, the coverage of the layer 
too + 1 fills exactly the vacancies of the layer mo — 1, and 
all layers k < m can be filled with the coverage of a 



corresponding layer with too < k < 2too. If mo is large 
enough compared to the roughness of the film, the layers 
with k > 2too have the coverage 0. The total specular 
intensity at the anti-Bragg point scattered by the rough 
film is equivalent to the intensity scattered by too — 1 lay- 
ers with the coverage 1 (corresponding to Eq. IA4|) , and 
by the layer too with the coverage 0.5. Therefore, A9 
always equals 0.5. 
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